Active gel model of amoeboid cell motility 
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We develop a model of amoeboid cell motility based on active gel theory. Modeling the motile 
apparatus of a eukaryotic cell as a confined layer of finite length of poroelastic active gel permeated by 
a solvent, we first show that, due to active stress and gel turnover, an initially static and homogeneous 
layer can undergo a contractile-type instability to a polarized moving state in which the rear is 
enriched in gel polymer. This agrees qualitatively with motile cells containing an actomyosin-rich 
uropod at their rear. We find that the gel layer settles into a steadily moving, inhomogeneous state 
at long times, sustained by a balance between contractility and filament turnover. In addition, our 
model predicts an optimal value of the gel-susbstrate adhesion leading to maximum layer speed, in 
agreement with cell motility assays. The model may be relevant to motility of cells translocating 
in complex, confining environments that can be mimicked experimentally by cell migration through 
microchannels. 
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I. INTRODUCTION 

Cell motility plays a role in key physiological processes such as wound healing, morphogenesis, and immunological 
response. In the past few decades it has been a main focus of cell biology, leading to the identification of many of the 
molecular players involved. In particular, the actin cytoskeleton, a network of polar, semiflexible protein filaments has 
been shown to be an essential part of the motility machinery. This network is highly dynamic and out-of-equilibrium: 
filament polymerization and depolymerization give rise to spontaneous network flows and molecular motors such as 
myosin-II interact with actin, behaving as active crosslinkers and exerting stresses on the network. 

The actin cytoskeleton has also attracted much attention in the physics community, since it represents an important 
example in the class of active living matter, that also includes, for example, bacterial suspensions, microtubule-kinesin 
solutions, and bird flocks [T]. A defining trait of these systems is a constant source of energy input, leading to non- 
equilibrium behavior such as pattern-forming instabilities and collective motion. The actin cytoskeleton, driven by 
hydrolysis of ATP, has been the focus of active gel theory [2H5] , a hydrodynamic approach providing a framework 
for quantitatively understanding cell biological phenomena such as the formation of contractile rings during cell 
division [5] , cortical flows during development [7] , and, of course, cell motility [5HTU] . 

Work on cell motility has focused on two basic migration modes. In the first mode, sometimes referred to as the 
mesenchymal mode, cell crawling on a flat surface is powered by the lamellipodium, a thin, fan-shaped structure in 
front of the cell body that is dense with actomyosin [11]. Crawling motility on 2D surfaces is understood as follows: 
polymerization of actin at the lamellipodium leading edge along with anchoring of new filaments to the substrate via 
focal adhesions generates a pushing force against the cell membrane, while contractility of the actomyosin network at 
the rear pulls the cell body forward. Steady moving states are maintained by turnover of actin filaments, a process 
known as treadmilling. Polarization of the actomyosin gel, strong adhesion, and high cell shape anisotropy of the 
lamellipodium are characteristic of this type of motility [T2] . Theoretical modeling based on one-component active gel 
theory has successfully accounted for several features of this migration mode, including the dependence of cytoskeletal 
flows on contractility and cell speed on filament turnover [8]. 

However, the environment that a cell usually encounters in vivo is a 3D extracellular matrix, providing only weak 
attachment points that the cell can push off of, thus making the lamellipodium mode not well suited. In fact, several 
recent studies have shown that cell types such as fibroblasts [T3] , leukocytes [H] , and cancerous human breast cells [TS] , 
can migrate in complex 3D geometries via an amoeboid mode. In addition, tumor cells can undergo a mesenchymal- 
to-amoeboid transition in the presence of protease inhibitors that prevent extracellular matrix degradation, thus 
favoring amoeboid- type movement which enables the cells to squeeze through small spaces [161 117j . In the amoeboid 
mode, there is often no leading edge actin polymerization, and motility is strongly dependent on myosin II-driven 
contractility of the actin network, generating cytosolic fluxes and often producing blebs. Despite growing interest in 
amoeboid motility, there has been far less theoretical investigation done into this mode than into its lamellipodial 
cousin. Contractility-induced permeation of the cytosol through the actin cytoskeleton is a critical aspect of this 
migration mode, and a theoretical description requires treating the cytoskeleton as a multi-component, poroelastic 
gel. This approach has successfully been used to quantitatively describe blebbing [18]. We propose here a model of 
amoeboid cell motility based on multi-component active gel theory [5]. 

Precisely how cell polarity and directional movement are maintained remain open questions. Polarity in 2D crawling 
cells is thought to be sustained by a complex signaling network involving localization of Rho family proteins to the 
lamellipodium leading edge [19 . In contrast, establishment of polarity in cells executing amoeboid motility is less well 
understood. For instance, polarity-associated signaling molecules are delocalized during 3D fibroblast migration [13] . 
Yet, a recent model of round breast tumor cells migrating in extracellular matrix |15j has uncovered a means for the 
cell to migrate persistently without a preexisting leading edge: myosin Il-induced contractility generates an instability 
of the homogeneous actin cortex of a stationary cell, breaking the spherical symmetry, and giving rise to cortical actin 
flows and cell motion |10j . This results confirms the critical role played by contractility in the spontaneous symmetry 
breaking of static round keratocytes [5D], keratocyte fragments [2T], and active droplets [52] . Nonetheless, it is not 
yet clear how the two fundamental processes involved in cell motility, namely, active contractility and turnover of the 
actin cytoskeleton, relate to the persistence and speed of cell motion. 

In this paper, we develop a model of amoeboid cell motility, treating the cell as a confined layer of a multi- 
component, poroelastic, active gel in a one dimensional geometry that is relevant to motility assays of cells confined 
in microchannels 0123]. We first elucidate the linear instabilities of a homogeneous, stationary gel, which results in 
a polarized, moving state. We find self-propelled moving steady-states at long times, and study the relation between 
motility, contractility, adhesion to the substrate, and filament turnover. 

Several recent theoretical studies have found pattern- forming instabilities in active gels of infinite extent [24H28] . By 
considering a finite film geometry, we see how these instabilities, coupled frictionally to a substrate, yield center of mass 
motion. Cell motility is a vastly intricate problem, involving feedback between biochemical signaling and mechanics. 
However, by following the modeling approach developed here, based on conservation laws and force balance, we are 
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able to uncover basic dependencies that necessarily remain in more complex models. 



II. MODEL OF AMOEBOID CELL MOTILITY 



Eukaryotic cell motility, whether lamellipodium-based crawling of keratocytes on a 2D surface [IT] or translocation 
of quasi-spherical cells through a 3D extracellular matrix [15] . depends on polymerization and depolymerization of 
the actin filaments and on active stresses within the actin network. Here, we model confined cell motility as follows. 
Consider a layer of an isotropic polymer network in a solvent held between two flat surfaces. The layer is initially 
homogeneous and at rest with respect to the surfaces. The length of the layer is L and its height is h. To represent 
loads, such as the nucleus trailing the motile apparatus or a viscous load encountered by the cell anterior, we introduce 
two identical "pistons" at each end of the gel layer; see Fig. [T] This geometry is in fact relevant to motility assays of 
cells confined in microchannels [5J [23] ■ 




FIG. 1. A layer of moving active gel of length L. The gel consists of polymer filaments (red lines) permeated by a solvent (light 
gray); it is bound above and below by two fixed horizontal surfaces, and on the left and right by two moveable pistons. The 
blue curve, determined by numerically solving Eq. |8| , indicates the polymer volume fraction in the steady-state of the layer 
moving towards the right, and illustrates that the layer is enriched in polymer at its rear. The parameter values used in the 
steady-state solution of Eq. Q are (fio = 0.5, £ = 33, kd = 200, a = 100, and = 0.46. The corresponding steady state speed 
of the film is Uoo = 16.7. 



A. Active gel description 



Active gel theory provides a continuum framework for studying the long time and length scale flows of the cy- 
toskeleton [2J 2] [5] . In this theory, the cytoskeleton is described as an actin polymer network permeated by a cytosol. 
The network is cross-linked by several types of proteins such as a-actinin and filamin; rheological studies have shown 
that actin networks undergo solid-to-liquid behavior on timescales greater than 1-10 s [25]. Moreover, actin filaments 
interact with molecular motors such as myosin II, that, via ATP hydrolysis, exert active stresses on the network. 



1. Polymer mass balance 



In our model, gel motion arises from active stress-induced gel phase separation and network turnover. To see how 
this comes about we start with a simplified description of active gel theory, in which spatial dependence occurs only 
along the ir-direction. First, the polymer volume fraction 4>(x, t), whose value is 4>o m the homogeneous state, satisfies 
the continuity equation 

d t 4>- Ud x 4> + d x J p = -kd(<t>- <Ao) • (1) 

In this equation, U(t) and J p (x, t) are, respectively, the layer velocity and polymer flux with respect to the fixed frame 
and fcrf is the bulk filament depolymerization rate. In our description, filaments undergo uniform bulk polymerization 
with the rate kd <fio, assuming that diffusion of free monomers is fast enough that we may consider their concentration 
to be fixed at its homogeneous, unperturbed value. Our approach is different from those that restrict filament 
polymerization to one or more boundaries, thereby breaking the symmetry of the problem by hand. We note that the 
coordinate x is measured with respect to a fixed point in the frame moving with velocity U. Incompressibility of the 
combined polymer plus solvent system implies that the length of the layer is constant. 
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2. Forces on the gel 

Flows involved in cell motility occur at low Reynolds numbers and thus the gel layer is always at mechanical 
equilibrium. The forces acting on the polymer component of the gel include a force due to stresses in the network; 
an osmotic force due to composition inhomogeneities; polymer-solvent friction; and polymer-substrate friction. The 
polymer-solvent friction force is much smaller than the polymer-substrate force: the ratio of the two is Tj s /(^£ 2 ), 
where r/ s is the solvent viscosity, £ is the polymer-substrate friction coefficient, and £ is the network mesh size. Taking 
T] s = 10~ 3 Pa.s, £ = 10 16 Pa.s/m 2 [8] (for a layer thickness of 1 /xm), and £ — 50 nm, we find rj s /(^£ 2 ) ~ 10 -4 , and 
thus we may neglect momentum exchange between polymer and solvent. As a result, the balance of forces on the 
polymer component of the gel gives 

d x (a n - n) = £J p , (2) 

where a D is the xx-component of the network stress and II is the osmotic pressure. We note that in Eq. ^ the 
polymer-substrate friction force is written as £J p = ^4>v pi as opposed to the usual form £v p , where v p is the polymer 
velocity, to account for the dependence of the friction force on polymer volume fraction. 

The network stress, tr n , contains passive and active contributions. Here, we consider the liquid limit of the gel, 
valid on timescales greater than the Maxwell time tm = 1 — 10 s. In this limit, the passive part of the network stress, 
arising from polymer convection and crosslink remodeling, has the viscous form rj p d x J p , where rj p is the polymer 
viscosity JS] 1301 131] • However, this contribution does not qualitatively affect the flows and motile behavior of the 
gel layer and we do not consider it further [32]. The active part, on the other hand, is essential for motility. Since 
the magnitude of the active stress, arising from motor activity on filaments, increases with <j) a simple choice for the 
network stress is [33] 

a n = C0- (3) 

A linear dependence of the active stress on polymer concentration has also been considered in recent work on motile 
active droplets [22]. We consider positive values of the activity coefficient (, corresponding to contractile behavior of 
active actin networks [3HI3S]- Moreover, £ > can give rise to contractile instabilities and gel phase separation [2"7ll34| . 
We also note that the active stress has been interpreted in the context of blebbing and amoeboid motility as a myosin- 
induced active hydrostatic pressure [TBI [3fi] . Though this active stress results from the isotropic part of the stress 
tensor 5J, it is physically meaningful since the polymer component of the gel is generally compressible. 

The osmotic pressure, II, is formally given by II = —F/V + (pSF/S(j), where F(<p,d x (j>) is the ^-dependent part 
of the gel free energy and V is the gel volume [3TJ- Physically, II acts to saturate the linear instability causing 
gel phase separation and to smooth the interface between polymer-rich and polymer-poor regions. Thus, a simple, 
phenomenological form for II is 

n = a(0-^ o ) 3 -7^0, (4) 

where a and 7 are positive coefficients [37]. A term linear in (f> in II, describing filament diffusion, has been omitted 
since it can be absorbed into er n . Finally, combining Eqs. Q-Q we obtain the equation of motion for <f> 

d t 4> - Ud x cf> + -dl{W - ol{cI> - O ) 3 + idlcfy = -k d {<j> - <f> ) . (5) 

Equation ^ has, minus the convective term —Ud x (p, the form of a Cahn-Hilliard equation including a chemical 
reaction term [38 . This model has, in particular, been used to study spinodal decomposition and steady states of 
binary solutions with chemical exchange between the components |39t 1 10|. By including the convective term, we will 
show how Eq. ^ predicts inhomogeneous, moving steady-states. 



3. Boundary conditions and global force balance 

The solution of Eq. ([5]) requires the specification of four integration constants, determined from boundary conditions 
at the layer ends, and the unknown layer velocity U (t) . No- flux boundary conditions on the polymer are expressed 
as J p (0,t) = Ucj)(0,t) and J p (L,t) = Ud>(L,t). These equations imply no-flux of the solvent at the boundaries, as 
well. In addition, assuming that there are no specific surface energies or dissipation at x — and x — L implies that 
d x (j)(0, t) = and d x (j)(L, t) = [U22]. Finally, U(t) can be found by requiring that the gel layer plus the end pistons 
are at mechanical equilibrium, namely 

&U + t [ J p dx = 0, (6) 
Jo 



5 



where is the friction coefficient between the piston loads and the surface. 



B. Contractile instability, phase separation, and moving states 



1. Linear instability 

In this section, we demonstrate that an instability of the homogeneous, stationary gel layer gives rise to gel phase 
separation and layer motion. We will focus on the onset of unstable states, depending on the activity parameter 
£, the depolymerization rate kd, and the drag coefficient £/. First, a straightforward linear stability analysis can be 
performed in the limit £; — > oo, i.e., for fixed pistons, in which case the eigenmodcs 4>(x, t) — cf>o ~ e A ™ * cos (nirx/L), 
n = 1, 2, 3, . . ., substituted into Eq. ([5]), yield the dispersion relation 

i/n7 ^ 2 (c-(^\y (7) 



£ \ L J \ \ L J 

Equation ^ reveals that activity drives the linear instability of the gel; the mixing effect of filament depolymerization 
suppresses long wavelength modes; and the smoothing term —jd 2 <p cuts off short wavelength modes. 

The stability analysis for the case of moving pistons is more involved, as the third spatial derivative of the eigenmodes 
no longer vanishes at the layer ends. To highlight the role of £, kd, and £j, we introduce the dimensionless quantities 
C = (L 2 /"/, k d = fc d ^i 4 /7, £ ; = £;/(££), a = aL 2 /j, J p = J p ^L 3 /j, and U — U^L 3 /^. Furthermore, we non- 
dimensionalize x by L and t by £L /7, keeping the old variable names for the new ones. As a result, the dimensionless 
equation of motion for <f> is 

d t <t> - Ud x + d 2 x (U - 5(0 - O ) 3 + d* x <P) = -kd{<t> - <Po) ■ (8) 
Writing <fi(x,t) — 0o — e* nt <f>n(x), inserting this ansatz into Eq. Q, and keeping only linear terms in <p n yields 

C» +<€(X) = -(k d + \n)<j>n(x) , (9) 

where primes (') denote differentiation with respect to x. Writing the layer velocity as U(t) — U n e Xnt , the boundary 
conditions on the third derivatives of <f> and the flux J p at x — and x = 1 read 

4>' n {o) = o, 0'„(i) = o; (ioa) 
^'(o) = u n <p , C(!) = ^o; (10b) 



while the global force balance is 



&Un + (<f>, 



(11) 





The modes with non-zero U n are antisymmetric about x = 1/2. We thus obtain 

<j) n {x) = A n+ sm[k n+ {x- 1/2)] + sin [fc„_ (x - 1/2)] , (12) 

where A n ± are constants and 




(13) 

The condition for non-trivial values of A n ± is obtained by setting the determinant of the two-by-two matrix obtained 



from the boundary and force balance conditions, Eqs. (10 )-(!!), to zero. This results in the characteristic equation 



fc n+ cos(fc„+/2) ( %^cos(fc„_/2) - 2(C 2 - fc J 2 l _)sin(fc n _/2) 




fc„_ cos (fcn_/2) I '^T^ cos ( fc »+/ 2 ) - 2 (C ~ *»+) sin (fc„+/2) 1=0, (14) 



G 



which can be solved numerically for the growth rate A„ for given £, kd, and £/. In the following, we will focus on values 
of £ for which n = 1 is the most unstable mode, corresponding, approximately, to ir 2 < £ < 9ir 2 . However, in principle 
higher, odd modes are excited for larger £, though the lowest odd mode is likely the most relevant for cell motility. We 
note that a finite drag coefficient £; affects the dispersion relation in a non-trivial way, as can be seen from Eq. (14 1. 
Qualitatively, we expect that the smaller the piston drag, the more unstable the gel layer is. Indeed, Fig. [2] reveals 
that the growth rate of the first mode, Ai, decreases as £ ; increases. Interesting, if, in the polymer force balance, the 
polymer-solvent drag were to dominate the polymer friction with the substrate, i.e., J p — <f>U ~ d x (a n — II), then A 
would no longer depend on £j, since no-flux at x = 0, 1 would mean that the eigenmodes for (j>(x,t) — cf>o are always 
of the form cos (mrx). 
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FIG. 2. Growth rate of n = 1 mode versus piston drag. Ai versus is shown for three different activities, as indicated. The 
other parameters are kd = 190 and <f>o = 0.5. 




FIG. 3. State diagram in the (kd,C>^i) s P ace - Moving states exist in the region bound by the surfaces Ai = and kd = 0. 
The initial polymer volume fraction is 4>o = 0.5. 
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2. Gel phase separation and moving steady-states 



In the regime Ai > 0, a contractile instability of the gel layer appears, which is triggered by active stress, but is 
counteracted by filament turnover. However, filament turnover is essential for moving steady-states of the gel to exist. 
Indeed, if kd = 0, polymer conservation implies that at steady-state J p (x) = U<j){x); substituting this into the global 
force balance, Eq. (J6|, and solving for U , it follows that the layer speed is zero. 

An inhomogeneous polymer density profile and filament turnover give rise to moving states in the regime Ai > 0. 
The state diagram in Fig. [3] indicates the region in the (£, kd, £;) parameter space in which such moving steady-states 
exist, bound by the surfaces kd = and Ai = 0. The moving states are characterized by an inhomogeneous profile 4>{x) 
and a non-zero steady state velocity U^, which are obtained by numerically solving the non- linear Eq. Q, completed 
by boundary conditions and global force balance. Figure [l] shows an example of a moving steady-state for non-zero 
kd that develops after the initial instability of the gel, and in which the layer is enriched in polymer at its rear. The 
polymer volume fraction, <j), is indicated by the blue curve in Fig. [I] For non-zero kd, net filament depolymerization 
occurs in parts of the layer where <j) > 4>o, whereas net polymerization occurs in regions where <f> < 4>q. As a result, 
polymer flux is directed from right to left in Fig.[T] exchange of momentum between the gel and the substrate provides 
a rightward force and drives motion of the layer to the right. 

The moving states can be characterized by the load force - vel ocity relationship for the moving layer, described here 
by the dependence of £ z on U^; see Fig. |4j Figures 4(a) and 



4(b) show that £ ; decreases with increasing U^ for, 



respectively, fixed C and varying kd and fixed kd and varying (. We note that if the dispersion relation is such that 
X 1 > (obtained for £ ; — > 00), then one must have — > 00 as Uoo — » 0, since only rigidly fixed pistons can prevent 
a linearly unstable gel from being set into motion. On the other hand, if A x < (obtained for £ ; — > 00), moving 
steady-state exist for < £ ; , where depends on £ and kd- In fact, £ ; is the load friction for which the slab is 
marginally stable: Ai(£, fcrf,^) = 0. The state diagram in Fig. [3] shows that if the depolymerization rate is kd is 
large enough, the mixing effect of filament turnover suppresses gel phase separation and self-sustained moving states 
no longer exist. This illustrates the subtle role played by filament turnover in our model: it is necessary for moving 
states to be sustained, yet if it is too high, cell polarity, needed for motion, is abolished. 




U 



FIG. 4. Load friction coefficient £ ; versus steady-state layer speed Uao. (a) £ ; vs U^ is shown for three indicated values of 
kd and C, = 30. (b) £ ; vs Uoo is shown for three indicated values of C, and kd = 190. £ ; tends to a finite value as U 00 — > if 
the linear growth rate Ai < 0. In (a) and (b), the parameter values used in obtaining the steady-state solution of Eq. (JsJ) are 
<j>o = 0.5 and a = 100. 



Motion of the gel layer requires a transfer of momentum from the substrate to the gel. The dependence of the 
steady-state layer speed, Uoo, on the polymer-substrate friction coefficient, £, allows a comparison between our model 
and cell motility experiments that probe the dependence of migration speed on surface adhesion. Figure [5] reveals that 
there is an optimum value of £, which maximizes Uoo- This prediction agrees with biphasic behaviour of migration 
speed as a function of surface adhesion observed in lamcllipodium-driven crawling cell motility on flat surfaces [431 144] 
and in amoeboid motility of cells confined in microchannels |23j . Our model predicts that the optimum value of £ 
shifts to higher values as the activity £ increases, in agreement with the work of Ref. [44]. In the model presented 
here, the layer speed is small for low £ because the polymers do not get enough traction from the substrate, and is 
small again for high £, since, according to Eq. ([5]), the amplitude of gel phase separation, and hence the layer speed, 
decreases with decreasing 
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FIG. 5. Steady-state layer speed as a function of polymer-substrate friction. The dimensionless speed Uoo/(kdL) versus 
dimensionless friction coefficient £ kdL 4 /^ is shown for three values of activity, as indicated. The other parameters are £j / (£L) = 
100, a = 100, and <f> = 0.5. 



We have developed a model of amoeboid cell motility that accounts for observed characteristics of this migration 
phenotype such as increasing cell speed with increasing contractility and a non-monotonic dependence of cell speed 
on surface adhesion. The simple approach used here puts into quantitative relief the important roles played by con- 
tractility and filament turnover, considered as the two driving agents of cell motility. While we found the dependence 
of migration speed on active stress to be consistent with experiments on 3D motile cells subject to varying motor 
activity, our predictions involving the filament turnover rate are more nuanced. Theoretical work has underlined the 
importance of the actin depolymerization rate on speed [H H5J , yet our work suggests that this rate may be tuned to 
maintain a balance between filament recycling, needed for motion, and cell polarity, involving actin density variations. 

Amoeboid cell motility involves cytosolic permeation through the actin cytoskeleton generated by molecular motor 
action. In addition, the polymer density in the gel is non-uniform, as the cell rear is generally enriched in actin 
filaments. A full accounting of these two features requires a multi-component active gel approach, which we have 
pursued here. We note that related, multi-component models have been developed in other contexts (IHIHJ- Our 
model, treating the cytoskeleton as an active isotropic gel undergoing delocalized filament turnover, is relevant to 
bulk and cortical cytoskeletal flows involved in amoeboid movement. There is, at this stage, no evidence of long range 
orientational order associated with amoeboid migration, which might be related to the absence of leading edge actin 
polymerization in this phenotype 48, 49J. This is very different from the situation encountered in lamellipodium- 
based cell motility, in which leading edge polymerization sets a preferred filament ordering direction, perpendicular 
to the leading edge |50| . The motility mechanism that we describe here also has similarities with the one proposed 
in |10j , which was based on contractile instabilities of the actomyosin cortex generating steady state cortical flows; our 
present work, in fact, extends these results to cell types where the actomyosin system is not localized on the membrane 
but distributed in the cell bulk. We believe that such mechanisms based on actomyosin contraction localized at the 
back of the cell and filament recycling could in fact be very general in the context of cell motility in confinement, 
which can be now routinely studied in vitro in microchannel assays [HI H3] • 

We point out finally that the modeling approach developed here may be extended to describe blebbing cell motility, 
a subclass of amoeboid motility involving coordinated blistering of the cell membrane and motion of the cell [36] - 
Cortical contraction and squeezing flow of the cytosol are implicated in this type of motion. An additional ingredient, 
which we have not considered here, is the convection of dissolved motors and actin monomers by the cytosol into the 
advancing bleb. Further work is needed to understand how blebbing, cytosolic flow, and protein transport conspire 
to produce cell motion. 
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